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The projection-operator formalism of Feshbach is apphed to resonance scattering in a single- 
channel case. The method is based on the division of the full function space into two segments, 
internal (localized) and external (infinitely extended). The spectroscopic information on the reso- 
nances is obtained from the non-Hermitian effective Hamilton operator HcS appearing in the internal 
part due to the coupling to the external part. As is well known, additional so-called cutoff poles 
of the S matrix appear, generally, due to the truncation of the potential. We study the question 
of spurious 5* matrix poles in the framework of the Feshbach formalism. The numerical analysis is 
performed for exactly solvable potentials with a finite number of resonance states. These potentials 
represent a generalization of Bargmann-type potentials to accept resonance states. Our calculations 
demonstrate that the poles of the S matrix obtained by using the Feshbach projection-operator 
formalism coincide with both the complex energies of the physical resonances and the cutoff poles 
of the 5* matrix. 
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I. INTRODUCTION 

The Feshbach projection operator (FPO) formal- 
ism [H, [l] is a powerful method for the description of 
resonant scattering and reactions involving light nu- 
clei In recent years the FPO technique has been 
applied to numerous other systems like quantum dots 
and microwave cavities [1, H, 0, H, H, Uo', Tl*! and atoms 
in a laser field [H, [13, [lj| • In its original formulation, 
the formalism is based on the introduction of projec- 
tion operators Q and P, QP = PQ = 0, P + Q = 1, 
which project, respectively, onto the discrete states 
of a closed system and the continuous spectrum of a 
reservoir when their interaction is neglected. Reso- 
nance states then naturally appear as bound states of 
the former closed system embedded into a continuum 
of open channels, due to the coupling really existing 
between the closed system and the reservoir. In other 
words: the starting point of the FPO formalism is the 
assumption that the scattering event is confined to a 
certain compact part of the available space ^Bl • This 
region constitutes the so-called interaction region and 
can be described by the Q subspace. Outside this 
region (in the P subspace) the interaction is absent 
so that the motion of scattering fragments depends 
(apart from the total energy E) only on their internal 
states. Each combination of internal states of all frag- 
ments is called a channel of reaction since it specifies 
a set of configurations (depending on E) in which the 
system can be found long before and long after the 
scattering takes places. 

The FPO formalism exploits the concept of an ef- 
fective Hamiltonian H(.g to describe the open system 
resulting from the interaction between the idealized 
closed system and the reservoir. The operator HcS 
is, naturally, non-Hermitian and depends explicitly 
on energy. Its complex eigenvalues z\ are energy de- 



pendent. The solutions of the fixed-point equations 
for the eigenvalues provide approximately both energy 
positions and inverse lifetimes (widths) of resonance 
states [3]. Using the FPO formalism, an expression for 
the S matrix can be derived ^ . It contains the com- 
plex eigenvalues z\ with their full energy dependence. 
The energy dependence is important especially in the 
neighborhood of decay thresholds and in the regime 
of overlapping resonances [H, [13] ■ 

As it is well known, analytic properties of the S ma- 
trix are very sensitive to both the detailed form of the 
potential and the behavior of the potential at infinity. 
Even at a point where the potential equals to zero 
with a computer precision, the truncation strongly 
affects the picture, especially the number of the S- 
matrix poles. In the simple case of the scattering by a 
potential of a compact support, i.e., a potential van- 
ishing outside a given cut off radius, the S matrix has 
an infinite number of discre te p oles in the lower part 
of the complex k plane [H, [l^ [lO] whereas an expo- 
nential asymptotic form of the potential can lead to 
a finite number of poles (see, e.g., [HI). This means 
that for a truncated potential one obtains not only 
physical (resonance) 5'-matrix poles but also so-called 
cutoff poles These poles are not an artefact; they 
are correct poles of the truncated potential but do not 
cause the characteristic phase shift by tt. Therefore 
there is a need to distinguish between the two kinds 
of poles [13] ■ One way to do so is to use the fact that 
the positions of the resonance poles are not affected 
by changing the model parameters such as, for in- 
stance, the cutoff radius. From a mathematical view- 
point both types of poles are correct but they have 
different origins. This is the reason why the problem 
of separation of cutoff poles from the resonance poles 
is widely discussed in the literature devoted to the S 
matrix. 
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To the best of the authors' knowledge the emergence 
of spurious complex eigenvalues has received little at- 
tention in the context of the FPO formalism (23l. [23|. 
We think that the reason for that may reside in the 
fact that some authors (e.g., [1^) find the concept of 
the non-Hcrmitian effective Hamiltonian unsuitable in 
the case of potential scattering. In our opinion, how- 
ever, the problem is not investigated in necessary de- 
tails and the current paper is just devoted to fill in 
this gap. 

To avoid unnecessary complications we apply, in 
the present paper, the FPO formalism to resonance 
scattering in its simplest form by considering only 
single- channel (elastic) s-wave scattering. In doing 
so we exclude the appearance of the Feshbach (core 
excited) resonance states that are present in the gen- 
eral case (see, e.g., [11|). Thus only potential (shape) 
resonances may appear in our approach. We assume 
the potential to have a finite number of resonance 
states and the continuous spectrum to fill the posi- 
tive semiaxis resulting as solutions of the usual radial 
Schrodinger equation. We construct exactly solvable 
potentials by applying the method of supersymmetric 
quantum mechanics (SUSY QM) |27| to the inverse 
scattering problem (see, e.g., Refs. (211 [28jV By ex- 
act solvability we mean the situation when solutions 
of the Schrodinger equation are available in an ex- 
plicit form and, in particular, are expressed in terms 
of elementary functions. 

It is the authors' opinion that the potentials, the 
scattering data of which we set ourselves, represent 
a good testing ground for numerous schemes of reso- 
nance calculations. Such calculations show the sub- 
stantial relevance of the concept of the non-Hermitian 
effective Hamiltonian to resonance scattering on a fi- 
nite range potential. We show that in this case the 
FPO formalism gives very accurate results both for 
the scattering phase shift and the positions and widths 
of physical resonances as well as for the cutoff poles 
of the scattering matrix. We discuss the fitness range 
of the fixed-point approximation and omit from our 
discussions the question of why the estimation of posi- 
tion and width of the resonance states by this method 
might be meaningful. In calculations for concrete re- 
actions by using the FPO method, this approximation 
is never used since the S matrix contains the energy 
dependent functions z\{E). It does not contain the 
energy-independent values that characterize the posi- 
tions and widths of the resonance states and are ob- 
tained from, e.g., the solutions of the fixed-point equa- 
tions [l6!,[i3|- present paper, we determine the 
poles of the S matrix exactly within the FPO formal- 
ism and compare them with the results of the exactly 
solvable potentials. 

The paper is organized as follows. In Sec. |TT]we re- 
call the methods of SUSY quantum mechanics in the 
context of scattering theory and construct Bargmann- 
type potentials supporting resonance states. Using 
the standard procedures of quantum scattering the- 



ory in Sec. IIIIl we calculate the S'-matrix poles for 
truncated Bargmann-type potentials. In Sec. llVi we 
write down the basic equations of FPO formalism used 
in the paper. Our numerical scheme for the computa- 
tion is mainly based on the consideration of quantum 
scattering on billiards with one attached lead @ . This 
model is formulated in the tight-binding approxima- 
tion (Anderson model). In Sec. |V]we present numeri- 
cal results obtained by using the two considered meth- 
ods. In the last section some conclusions are drawn. 



II. BARGMANN-TYPE POTENTIALS 
SUPPORTING RESONANCE STATES 

As it is known, the SUSY approach, when re- 
stricted to the derivation of new exactly solvable quan- 
tum problems, is basically equivalent to the Dar- 
boux transformation method (see, e.g., Ref. [29j). 
Therefore we use SUSY and Darboux transformations 
as synonyms. The whole class of potentials known 
as Bargmann-type potentials (see [30|) may be ob- 
tained from the zero potential with the help of ei- 
ther usual SUSY transformations or their confluent 
forms (see ,31, 32, 33]). Typically, such a potential 
has an exponentially decreasing tail and supports a 
finite number of bound states and no resonances. Its 
S matrix (as well as the Jost function) is a rational 
function of the momentum k (see, e.g., |20|). Any such 
potential may be obtained by a proper chain of SUSY 
transformations with real factorization constants (see 
the next section for details). Nevertheless, as it is 
shown in ^34] , the use of complex factorization con- 
stants in higher order transformations [1^ (see also 
(36j ) may lead to a real potential and corresponds 
to an irreducible supersymmetry. The use of such 
SUSY transformations permits us to enlarge the class 
of standard Bargmann potentials by including poten- 
tials supporting resonance states. 



A. Darboux transformation method 

In this section we shortly recall the definition 
and main properties of the Darboux transformations 
method necessary for subsequent analysis. The inter- 
ested reader can find a more detailed exposition else- 
where m, [m, m m H US n, [13, m, m, no, sn. 

In its pragmatic formulation the method essentially 
consists in getting solution ip of the (transformed) dif- 
ferential (Schrodinger) equation 

h^^ = E^, h, = -f^ + Viir) , (1) 
ar^ 

by applying differential transformation operator L of 
the form 

L = -d/dr + w (r) , (2) 
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to a known solution ijj of another (initial) equation 

hoi^^E^, ho + Voir), (3) 

If = Lip, corresponding to the same value of the pa- 
rameter E. Here the real- valued function w (r) called 
the superpotential is defined as the logarithmic deriva- 
tive of a known solution to Eq. ([3]) denoted by u 

w ^ u' (r) /u (r) , hou — au (4) 

with a < Eq, where Eq is the ground state energy 
of the Hamiltonian ho if it has a discrete spectrum 
or the lower bound of the continuous spectrum other- 
wise. Function u is called the transformation or fac- 
torization function and a is known as the factorization 
constant or factorization energy. The potential Vi is 
defined in terms of the superpotential w as 

(r) = Vo (r) - 2w' (r) . (5) 

The knowledge of all two-dimensional solution space 
of the initial equation with a given value of i? G C pro- 
vides the knowledge of all solutions of the transformed 
equation corresponding to the same value of E. In 
particular, if all solutions for E belonging to the spec- 
trum of ho are known (so-called physical eigenfunc- 
tions of ho) the method provides us with all physical 
eigenf unctions oi hi. 

Since the above procedure does not depend on a 
particular choice of the potential Vq, the transformed 
Hamiltonian hi can play the role of the initial Hamil- 
tonian for the next transformation step. In this way 
one gets a chain of exactly solvable Hamiltonians ho, 
hi, hn with the potentials Vo, Vi, . . . , Vn- To 
avoid any confusion we mention that everywhere, ex- 
cept for especially mentioned cases, we shall use sub- 
scripts to distinguish between quantities related to 
different Hamiltonians, ho, hi , ... and shall omit 
them when discussing general properties regarding all 
Hamiltonians. 

First-order Darboux transformation opera- 
tor Ljj+i, as defined by Eq. relates solutions of 
two Hamiltonians hj and /ij+i. If one is not interested 
in the intermediate Hamiltonians hi, hn~i and 
all factorization energies are chosen to be different 
from each other, the whole chain may be replaced 
by a single transformation given by an nth-order 
transformation operator, denoted by L'-"\ defined 
as a superposition of n first-order transformation 
operators. A compact representation of this operator 
is given by [l^j 

iJn(r,k) = L^^Uo{r,k) (6) 

= W (ui, ...,Un,lpO {r, k)) [Ui, ...,Un), 

where ipvi {r, k) is a solution to Eq. ([3]) corresponding 
to the energy E = k^ and ipn {r, k) satisfies 

hnipn {r,k)^ EiPn {r,k), E^k". (7) 



The transformation functions Uj , although labeled by 
a subscript, are eigenfunctions of the initial Hamilto- 
nian 

hoUj (r) — a'juj (r) . (8) 

These should be chosen in a way that the Wron- 
skian W{ui,...,Un) is nodeless and either real or 
purely imaginary for r g (0, oo). These conditions 
guarantee the absence of singularities in the potential 

Vn = Vo-2^lnWiui,...,Un) (9) 

defining the Hamiltonian hn of Eq. ([7]) and its real 
character for r S (0,oo). In particular, factorization 
constants should be either real or come in complex 
conjugate pairs with corresponding factorization solu- 
tions being either real or in pairs complex conjugate 
to each other. Formula ([6]) is valid for any E = k^ 
except for k = aj (j = 1, . . . , n). For these values of k 
the corresponding solutions are 

ipn {r,aj) = W'^^'> (ui, ...,Un) [ui, ...,Un) , 

(10) 

where W'-^^ {ui,. . . , Un) is the (n — l)st order Wron- 
skian constructed from ui , . . . , u„ except for Uj , 
j = l,...,n. 

B. Jost function for a special chain of 
transformations 

Let us choose the following set '4l'] of eigenfunctions 
of the Hamiltonian ho © as transformation functions 
for the Darboux transformation of order 2n: 

ui (r) , vi (r) , U2 (r) ,V2 (r) , . . . ,Un (r) , Vn (r) , (11) 

hnuj (r) = a]u^ (r) , hovj (r) = P^Vj (r) . (12) 

Here a's and /3's should be different from each other 
and chosen in a way that the corresponding factor- 
ization energies, if they are real, are smaller than the 
ground state energy of ho when it has a discrete spec- 
trum, or less than the lower bound of the continuous 
spectrum otherwise. No such restrictions are imposed 
on a's and /3's for complex factorization energies. 

We distinguish between the functions u and v from 
their behavior at the origin. The functions v are regu- 
lar, Vj (0) = 0, and hence are uniquely defined up to a 
constant factor, that is not essential for our purpose. 
The functions Uj are irregular at the origin, Uj (0) ^ 0, 
and form the so-called singular family. 

In [41] it was shown that the chain of transforma- 
tions with transformation functions (jlip transforms 
the initial Jost function Fo (fc) of Hamiltonian ho to 
the Jost function, 

Fn{k)^Fo{k)l[^-^, (3,^^b,, (13) 
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J with flj > 



No 



corresponding to the Hamiltonian /i„. Since the Jost 
function should be analytic in the upper half of the 
complex k plane (see Refs. [H M, llj, M M, S M, 
[2^) all &'s must be real and positive. This avoids the 
appearance of the so-called redundant poles, which 
occur as poles of the Jost fmiction or zeros of the S 
matrix. Every purely imaginary aj — ia 
corresponds to a discrete level Ej = —aj of h 
restriction, except the ones discussed above, is im- 
posed on a's. Complex a's coming in pairs with real 
parts of opposite signs may correspond to a visible 
resonance. Thus, we say that every pair of complex 
numbers aj — ±Re [ctj] -|-ilm [aj] corresponds to a res- 
onance and to mirror resonance states with complex 
energy, 

-Bf ' = (±Re [aj] + ilm \aj]f = (14) 
(Re [ujf - Im [ujf^ ± 2iRe [a^] Im [uj] . (15) 

In this case, in accordance with the analytic proper- 
ties of the Jost function (fT3| , the technique developed 
in i^lj] remains stable for Re [uj] < 0, and Im [aj] < 0. 
States with |Re[aj]| > |Im[Q;j]| correspond to vis- 
ible resonances provided Re [aj] — Im [aj] is small 
enough [isj . 

In the one-channel case the S matrix is a single- 
valued function of wave number k. 



S 



F{k) 



(16) 



In our approach (fc) differs from Fq (k) by a ratio- 
nal function of momentum k. Therefore the expres- 
sion for phase shift S„ (fc) becomes rather complicated 
when the number of transformation functions is suffi- 
ciently large. An alternative expression, which is more 
convenient for practical calculations, is ^l| 



— arctan I — 



Sn — ^0 — arctan 

(17) 

In Sect. |V] we apply the above described technique 
to obtain potentials with either one or two resonance 
states. Moreover, since our starting potential equals 
zero, when the usual technique of Darboux transfor- 
mations [41^ reproduces Bargmann potentials [soi] , the 
potentials from Sec. |V] are their generalizations to 
describe resonance scattering. Below we will denote 
them by Vbtq- It is worthwhile to note that the solu- 
tions of the Schrodinger equation for these potentials 
are expressed in terms of elementary (trigonometric) 
functions. 



boux method we investigate the scattering on trun- 
cated Bargmann potentials Vcut, i-e. on potentials 
equal to Bargmann potentials Vsrg for r < Rcut, and 
equal to zero for r > Rent- In this case, as it is well 
known, the analytic continuation of S (fc) to the com- 
plex k plane is a meromorphic function with infinitely 
many poles. In order to calculate their positions and 
widths we use the methods of usual quantum scatter- 
ing theory (see, e.g., [20l.f44|). 

For r > Rcut, where Vcut (f) — 0, the solution of 
Schrodinger equation ([T]) is a linear combination of 
plane waves, which we write as 

*r>i?e„t = cos (Scut) sin (fcr) -I- sin (Scut) cos (fcr) . 

(18) 

We denote the phase shift for the truncated potential 
as Scut ■ The phase shift Scut is obtained by solving Eq. 
([T|) for function ^'r<fl^„t (r) in the region r < Rcut and 
matching it to have form ifTS]) at r = Rcut- Function 
4'r<_R,<.„t(r) subject to the Dirichlet boundary condi- 
tion 5'r<i?cut (0) = 0, is uniquely defined (up to an 
inessential constant factor). Although at r ~ Rcut 
both '^r<Raut ^-nd (M! r<Rcut/ dr must be continuous, 
it is sufficient for our purposes to impose continuity 
on the logarithmic derivative (see, e.g., [3|), 





1 








J dr 




1 








t dr 



r=R^. 



r=R^. 



(19) 



which is independent of a multiplicative constant. 
From here we find the phase shift of the cutoff po- 
tential 



tan [Scut (k)] 



k cos {kRcut) - 7 sin {kRcut) 



k sin (kRcut) + 7 cos (kRcut) 
which, upon using Eq. (I16p . gives its S matrix 

k — ij 



Scut (k) 



t(fe) ^ g-2ifci?,„t 



(20) 



(21) 



The poles of Scut (k) are the roots of the transcenden- 
tal equation 



J — ik 

which we solve numerically in Sec. |Vl 



(22) 



IV. FESHBACH PROJECTION OPERATOR 
APPROACH TO POTENTIAL SCATTERING 

A. Basic relations of FPO formalism 



III. S'-MATRIX POLES FOR TRUNCATED 
BARGMANN POTENTIALS 

In order to compare the results obtained by means 
of the FPO technique and those obtained by the Dar- 



As was mentioned in the Introduction, in the FPO 
formalism [l|, Q the full function space is divided into 
two subspaces: the Q subspace contains all wave func- 
tions that are localized inside the idealized closed sys- 
tem and vanish outside of it while the wave functions 
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of the P subspace are extended up to infinity and 
vanish inside the system; see [l^ [l3 • This division is 
carried out by using the projection operators Q and P 
{QP = = PQ,P + Q = 1). The wave functions of 
the two subspaces can be obtained by standard meth- 
ods: the Q subspace is described by eigenf unctions of 
Hermitian Hamiltonian Ht that characterizes the lo- 
cahzed closed system with a discrete spectrum, while 
the P subspace is described by the states of Hermitian 
Hamiltonian Hf,, which has a continuous spectrum. In 
the FPO formalism, the closed system becomes open 
because of a really existing coupling between the lo- 
calized closed system and the reservoir, i.e., because of 
the coupling between the Q and P subspaces. Due to 
this coupling, some discrete states of the closed system 
become resonance states of the open system which, in 
general, have finite life times. 

In the present paper we are interested, above all, in 
the properties of the effective non-Hermitian Hamilto- 
nian of the open quantum system, which acts on the Q 
subspace and carries the influence of the P subspace. 
It reads 

Hcff = Hb + J2 ^bc^+^^Vcb . (23) 

c 

Here, _E+ — E + ie with e 0. Further, Vbc and Vet 
stand for the coupling operators between the Q sub- 
space (described by Ht) and the P subspace (envi- 
ronment, described by He). The operator Hcs is non- 
Hermitian and describes the localized system under 
the influence of the reservoir 

The non-Hermitian operator Hcs is complex- 
symmetric and depends explicitly on energy. Its eigen- 
values z\ and eigenfunctions 

(i?eff - za) = , (24) 

are complex. The eigenvalues provide not only the 
energies of the resonance states but also their widths. 
The eigenfunctions are biorthogonal. For more details 
see [l7| . 

We underline here that values like the S matrix and 
the cross section are independent of the manner how 
the Q and P subspaces are defined. However, in or- 
der to obtain the positions Ex and widths T\ of the 
resonance states from the eigenvalues z\ of HcS, the 
two subspaces have to be defined properly. Otherwise, 
the z\ have nothing in common with the spectroscopic 
values E\ — i/2 Fx oi the resonance states. This can 
be seen, e.g., from the fact that Fa -^OifQ + P^Q. 

Using the FPO formalism, the scattering matrix S 
can be written in terms of the effective Hamiltonian 
and the external scattering states \E, c) defined by 
{He - E) \E, c) = 0. It reads 0, [l^ 

5cc' = 5,,, - 2ni {E, c\ Vcb ^ Vbc' \E, c') . (25) 

Characteristic of Eq. ([^5]) is that it contains HcS 
with its explicit energy dependence. The energy de- 



pendence of -Heff plays an important role near de- 
cay thresholds and in the regime of overlapping reso- 
nances. The S matrix (|25p is always unitary. 

The FPO formalism may formally be considered as 
a certain generalization of the R matrix approach [45| . 
In both cases, the wave functions of the system are lo- 
calized in coordinate space {Q subspace in the FPO 
formalism) and coupled to an extended continuum of 
scattering wave functions (P subspace in the FPO for- 
malism). However, the standard spectroscopic param- 
eters of the R matrix approach do not contain any 
feedback from the continuum of scattering wave func- 
tions. In the FPO formalism, they are replaced by the 
energy-dependent functions E\ and Fa in which the 
feedback from the continuum is involved. 

We examine the concept of the effective Hamilto- 
nian in connection with the potential scattering on 
spherically symmetric potentials. In order to define 
the Q subspace that contains the localized part of the 
problem, we truncate the potential at a certain cutoff 
radius Rcut- The P subspace is defined then by the 
remaining part of the function space being extended 
up to infinity. The operators Vbc, Vcb describe the 
coupling between the two subspaces. In this paper, 
we consider a one-dimensional (ID) quantum system 
to which one lead is attached at a certain point. We 
will describe such a system in the framework of the 
tight-binding approach. 



B. One-dimensional tight-binding model for 
resonance scattering 

Let us consider the resonance scattering on trun- 
cated Bargmann potentials Vcut as described in Sees. 
[HI and ITTTl following the FPO technique. We choose a 
radius R such that the functions defined at < r < i? 
belong to the Q subspace, while the functions defined 
at r > P belong to the P subspace. In order to de- 
scribe the continuum (P subspace) properly, it should 
be P > Rcut- 

A common approach to solve the Schrodinger equa- 
tion in the context of the FPO formalism consists 
of the discretization of the spatial coordinate. The 
resulting matrix Hamiltonian is the so-called tight- 
binding Hamiltonian (see, e.g., 0]) which is widely 
used to model electronic transfer in molecules and con- 
densed matter. To obtain the matrix representation 
for the effective Hamiltonian HcS (l23|) . we choose a 
discrete lattice whose points are located at r = = 
ia, i = 1, . . . , iV (rAT = P) and approximate the sec- 
ond order derivative by the finite differences 

V'"«4{V'»-i-2V'» + ^.+i}, (26) 

where a = rj_i — is a lattice constant (being inde- 
pendent of i). Thus, for r = Ti, I = 1, . . . , A'' — 1, we 
obtain the following finite-difference (or tight-binding) 
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Schrodinger equation 

t{-tpi-i + 2ipi - V'i+i} + Uiipi = Etpi , 



(27) 



where Ui — Vcut{ri) and t = is the tight-binding 
couphng constant. 

Effects of scattering are introduced through the 
couphng between the box {Q subspace) and a semi- 
infinite lead (P subspace) that is attached at the 
point R = tm- In order to describe these effects, we 
present the sohition of the Schrodinger equation in the 
lead as 



-ikvi 



S'fpo (fc) e 



ikri 



i> N. 



(28) 



This equation defines the one-channel scattering ma- 
trix S'fpo (k)- In the above equation we use the stan- 
dard dispersion relation of the tight-binding model, 



where a is the lattice constant (see above) and E is 
the real energy of the system. At r = rjv, where the Q 
subsystem (box) is coupled to the P subsystem (semi- 
infinite lead), the Schrodinger equation takes the 
form 0, 113 



t {-ipN-i + (2 - e''='^)V'Ar} = E'iJjN-2ite-''''''" sm{ka) 

(30) 

Denoting ^/ = (Vji, i/'2, • ■ • , i/'w-i; V'w) , we get the 
matrix equation 



(S-Heff)*= (fffc-VF)] * = (31) 



£; = t [2 - 2 cos(A:a)] , 



(29) 



with coupling matrix Wij = Si^Sj^te'^''"' and 



f Ui + 2t -t 
U2 + 2t -t 
-< C/3 + 2t 









-t Un-1 + 2t 

~t 



-te 



ika 



2t / 



\ 2ite 



■ikrN , 



(32) 



i{ka) / 



Equation p2|) gives us the desired matrix repre- 
sentation for the effective non-Hermitian Hamilto- 
nian HcS- The Hamiltonian HoS obtained from the 
above pictorial derivation is completely equivalent to 
the overall Green function derived in [4^. The ma- 
trix equation PT|) describes the scattering on the dis- 
cretizcd Id quantum system with Bargmann potential 
at ri < Rcut and zero potential at > Rcut and with 
R = Rcut = rN ■ 

Although the poles of the S matrix are of no rele- 
vance for the scattering and reaction processes in the 
FPO formalism [l^ it is interesting to estimate 
their value. This can be done by using the following 
standard method First, the fixed-point equations 
for the positions of the resonance states in energy are 
solved, 

i?A = Rc[za]Is=b.> (33) 
and then the widths are defined by 

rA-2Im[zA]li;=£;, . (34) 

The solutions of these equations provide approxi- 
mately the energies E\ and widths Fx of the resonance 
states as long as the are small. In the following, 
we call the solutions of Egs. ([55|) and shortened 
the fixed-point solutions. 



In order to make a meaningful comparison of the 
results of the FPO method with those of the exactly 
solvable potentials, the poles of the S matrix obtained 
by using the FPO formalism should be determined 
exactly. This can be done by solving the equation 

Det [E - Heff (E)] = (35) 

in the complex E plane, which follows from the ex- 
pression l|25[) . provided that Vbc and Vet do not have 
poles. Equation psp differs from an eigenvalue equa- 
tion since the effective Hamiltonian depends on en- 
ergy. We found the solutions of Eq. (P5|) by deter- 
mining the intersection of contour lines for the zero 
values of the real and imaginary parts of the determi- 
nant. In order to distinguish the solutions of Eq. ([55)) 
from the fixed-point solutions, we call them S matrix 
poles in the following. 

Using the S matrix obtained from Eqs. (|3ip and 
(PSl) we calculate the phase shift according to its def- 
inition (|16p. 

Sppoik) ^ --In S'fpo (k) ■ (36) 

In all our numerical calculations presented in the 
next section, the lattice constant a is chosen to be 
0.01. We performed calculations also with a = 0.001. 
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These calculations require much more computer re- 
sources and change the results only slightly. The fixed- 
point equations are solved for R — Rent- 

V. RESULTS 
A. One-resonance potential 

To illustrate the differences between physical (res- 
onance) poles and unphysical (cutoff) poles more 
specifically we first apply the method described in Sec. 
im to construct a potential with one resonance state. 
This allows us to carry out a simultaneous compara- 
tive analysis of the smooth (nontruncated) potential, 
its truncated version (both obtained with the help of 
SUSY technique) and their cutoff counterpart result- 
ing from the FPO formalism. 

Our choice Vq {f) = in Eq. ([9]), which allows us to 
use solutions of the free particle Schrodinger equation, 
simplifies the calculations considerably. A resonance 
is obtained when two irregular transformation func- 
tions of type (|lip with real parameters ai and a2 are 
used in a chain of transformations, i.e., 

ui = exp (— iair) = cxp ((ai + m2) r) , 

U2 — exp (— ia2?') = exp ((ai — m2) r) . (37) 

As regular solutions we choose, in the current 
case, hyperbolic sine functions 

vi = sinh (6ir) , V2 = sinh (62'') • (38) 

The real constants a^, bi should be such that 6^ > 0, 
02 < ai < 0. The desired potential follows from Eq. 



([9]) where we have to calculate a fourth order Wron- 
skian W(ui,U2, vi,V2)- The calculations are simplified 
if we notice that this fourth order Darboux transfor- 
mation is equivalent to a set of two first order trans- 
formations and one second order transformation. If 
we choose for the first transformation the transforma- 
tion function u = ui, it results, according to Eqs. ^ 
and (O, in the zero potential difference (i.e., the ini- 
tial potential remains the potential of the free parti- 
cle). Then we have to change only the form of the 
three other transformation functions. These func- 
tions become L1U2, LiVi, L1V2 with L = Li given 
by Eq. ([2]) where w — —iai. Evidently, the func- 
tion L1U2 is, up to an inessential constant factor, 
the same exponential as U2, and the functions LiVi, 
L1V2 become proportional to hyperbolic cosines with 
shifted arguments. For the second transformation 
we choose the transformation function L1U2 which 
still does not add anything to the zero potential but 
changes the transformation functions LiVi —^ L2L1V1 
and L1V2 L2L1V2 [here L2 is given by the same 
formula ^ with w — —ia2], producing additional 
shifts in their arguments. After that we realize the 
second order transformation with the transformation 
functions 

L2L1V1 ~ sinh {bir - Ci) , L2L1V2 ^ sinh (62?" - C2) , 

(39) 

where 

tanhC» = ,2 2 ' ^ = 1'2 (40) 

which gives the desired one-resonance potential 



2{bl 



6| sinh (6ir — ^1) — 6^ sinh (&2»' — C2)' 



[62 sinh {bir ~ Ci) cosh {b2r - C2) - bi sinh (62^ - C2) cosh {bir - Ci)]" 

I 



(41) 



The potential Vires (see Fig. [U dashed line) rep- 
resents a generalization of a two-soliton potential de- 
fined on the positive semiaxis ^] . Instead of two dis- 
crete levels, present in the two-soliton potential, the 
potential ([4T|) has one resonance state. 

The Jost function reads 



Fires 



{k - ai){k - a2) 
[k + ibi){k + ib2) 



Hence the resonance occurs at = ai = —02 
with the mirror pole at k = a2 — 0-2 + iai, Ui 
i = l,2. 



(42) 



f tai 
< 0, 



The phase shift fTT]) is now given by 



Sir 



arctan ■ 



2aik 



arctan 



fc(6i 



(43) 

To be able to compare the results obtained by the 
two different methods, we draw all figures in the E 
plane [E = k^ for continuum and Eq. ([29]) in the 
tight-binding case]. 

First we compare the scattering phase shift cal- 
culated by means of FPO technique with that calcu- 
lated by Eq. (gS]); Fig. [2] The agreement of Eq. 
([55]) with the exact result (^5]) holds true up to high 
energy values and in a large range of cutoff radii. 
Even in the middle of tight-binding conductance band 
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separates from the chain of cutoff poles. It is chosen 
to lie close to the real energy axis. 



FIG. 1: One-resonance Bargmann-type potential 
(dashed line) at ai = —0.1, a2 = —2, fei — 1, and 
&2 = 2. Two-resonances potential (solid line) at 
ai = -0.1, 02 = -2, as = -0.08, 04 = -3, 61 = 0.2, 
62 = 0.1, 63 = 0.08, and 64 = 0.05. 



-1.5, 




4 6 
E 



10 



FIG. 2: Phase shift for one-resonance potential at 
ai = —0.1, 02 = —2, 61 = 1, 62 = 2, and i?cut = 5. 
Inset shows the difference between Sires and (5fpo. 



E — 2/t, where we cannot expect good agreement be- 
tween continuous and tight-binding models, the differ- 
ence i^ires— <^FPO does not exceed 0.1. The inset of Fig. 
[5] shows the difference between the values of the exact 
phase shift S and the numerical ^fpo- The maximum 
difference occurs near the resonance position. 

According to Eq. the phase shift is a sum 

of two terms. When the parameter values are those 
used in Fig. [2l both terms contribute with compara- 
ble weight even in the very neighborhood of the reso- 
nance. For other parameter values (e.g., ai — —0.01, 
02 = —2 for the resonance term and bi = 100, 62 = 
200 for the background), one term dominates in the 
neighborhood of the resonance and the phase shift is 
the standard one (i.e., tt) in this energy region. 

Let us now analyze the calculated spectroscopic 
data. As shown in [43] , a cutoff potential produces a 
chain of poles of the S matrix and there are no poles 
lying below this chain in the complex plane. A nar- 
row physical resonance of the nontruncated potential 
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FIG. 3: Poles of the truncated one-resonance poten- 
tial at fli = —0.1, 02 = —2, foi = 1, 62 = 2 and cutoff 
radius Rcut = 5. The symbols correspond to the roots 
of transcendental equation (|22p (o), the poles of the 
S matrix obtained from Eq. (I35|) (x), and the fixed- 
point approximation (|33}, (|34p {+). 

Figure [3] shows all poles (Re[i?] < 100) generated by 
the potential (|4ip truncated at Rcut — 5. The roots 
of transcendental equation ([22| are shown as circles, 
the crosses present the distribution of complex solu- 
tions of Eq. p5|) . and daggers stand for the results 
of the fixed-point approximation ([33]) . (|34|) . The reso- 
nance which is due to the nontruncated potential ()4ip 
is clearly separated from the cutoff poles. Its exact 
value is E = (-02 + iaif = 3.99 - iOA. The fixed- 
point approximation (|34p gives the width of the poles 
with a significant error. The S matrix pole (j35p re- 
produces this value to a high precision (relative error 
is less than 1% for shown poles). 

The picture related to the cutoff poles is, according 
to Fig. [31 the following: the S matrix poles ([55)1 co- 
incide with the roots of transcendental equation , 
whereas the widths (|34| obtained by using the fixed- 
point approximation arc far too small. That means 
that all poles of the S matrix are very good reproduced 
by the solutions of Eq. ([55]) . The fixed-point approx- 
imation (j33p allows one to accurately determine only 
the position of narrow resonances. This result agrees 
with the definitions (|33p and ([M)) according to which 
the fixed-point equation is solved only for the real en- 
ergy. 

In order to see how the resonance of the nontrun- 
cated potential separates from the chain of cutoff poles 
one can consider the dependence of the pole location 
on the cutoff radius Rcut- Indeed, only the physical 
resonance is almost independent of Rcut- AH the other 
poles move if the cutoff radius is changed. 

To show this dependence in detail we present the 
results from a set of calculations by using the FPO 
technique for different cutoff radii Rcut - First we con- 
sider the fixed-point approximation (Figs. [lland[n|). 
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FIG. 4: Solutions of fixed-point equations (|33|l and 
H34[l for one-resonance potential at ai = —O.I, 02 — 
—2, &i = I, and 62 ~ 2. Cutoff radius Rcut changes 
from 0.5 to 7 with step O.OI. With increasing Rcut, 
the trajectories move to small Ke{z), |Im(2;)| (with 
the exception of the spiralling trajectory). 

Or 
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FIG. 6: Trajectories of ^-matrix poles (|35p for one- 
resonance potential at ai — —O.I, 02 = —2, bi — I, 
and 62 = 2. Cutoff radius Rcut changes from 0.5 to 7 
with step O.OI. With increasing Rcut, the trajectories 
move to small Re(z), |Im(2:)| (with the exception of 
the spiralling trajectory). 



trajectory. 

The results shown in Figs. |4][6] demonstrate that 
the cutoff trajectories of the fixed-point solutions 

and of the S matrix poles ([55)1 depend strongly 
on Rcut- The spiralhng trajectories of the physical 
resonances arise from a lot of avoided crossings of 
the trajectories of neighboring cutoff trajectories. On 
this account, the spectroscopic values of the physi- 
cal resonances are influenced only a little by varying 
Rcut- The resonance location in Fig. [S]is stable with 
Rcut ^ 7 in contrast to those in Figs. [5] and [51 



FIG. 5: Solutions of fixed-point equations (|33p and 
l|34[) for one-resonance potential at ai = —0.2, 02 — 
—2, &i = I, and 62 ~ 2. Cutoff radius Rcut changes 
from 0.5 to 7 with step O.OI. With increasing Rcut, 
the trajectories move to small Re(z), |Im(z)| (with 
the exception of the spiralling trajectory). 



The trajectory of one of the eigenvalues has a strongly 
pronounced bight. The eigenvalue trajectory is spi- 
ralling around the correct value of the physical reso- 
nance {E — 3.99 — zO.4), see Fig. For the broader 
resonance {E = 3.96 — iO.8) shown in Fig. [5] the spi- 
ralling trajectory has less rotations and the localiza- 
tion of the resonance value becomes more difficult. As 
a result, the fixed-point approximation correctly indi- 
cates the position of the narrow resonance only. 

Figure [S] shows the solutions of Eq. ([55|l . They 
show a similar dependence on Rcut as the fixed-point 
solutions (Figs. |4] and [S]). However, the S matrix 
poles obtained by solving both Eqs. ([55)1 and 
coincide. The trajectory of one of the poles converges 
very quickly to the resonance value along the spiralling 



B. Two-resonance potential 

Let us now consider the more complicated case of 
a potential supporting two resonance states. To con- 
struct the Bargmann-type potential (see Sec. |TT| we 
are using eight transformation functions. As four ir- 
regular solutions from the set pT|) we choose two func- 
tions (|37p with 02 < ai < and the following two 
functions: 

U3 = exp {-ia^r) = exp [(03 + ia^) r] , 

U4 — exp {—ia^r) = exp [(03 — ia/C) r] (44) 

with 03 < a4 < 0. As four regular solutions we choose 
two functions (|38p and 

U3 = sinh (b^r) , V4 = sinh (&4r) (45) 

with all bi > 0. The desired potential follows from 
Eq. ^ where this time an eighth order Wronskian 
should be calculated. To simplify calculations we re- 
place, similar to the previous section, the eighth order 
Darboux transformation by a chain of four first or- 
der transformations, involving exponential transfor- 
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FIG. 7: Poles of the truncated two-resonance poten- 
tial at ai = —0.1, 02 = —2, 03 — —0.08, 04 = —3, 
bi = 0.2, 62 = 0.1, 63 = 0.08, &4 = 0.05, and cut- 
off radius Rcut = 5. The symbols correspond to the 
roots of transcendental equation (I22|) (o), the poles 
of the S matrix (|35|l obtained from the eigenvalues 
of i?eff (x), and the fixed-point approximation 

m (+)■ 
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FIG. 8: Solutions of fixed-point equations (I33|l and 
(|34p for two-resonance potential ai = —0.1, 02 — —2, 
03 = -0.08, 04 = -3, bi = 0.2, fe2 = 0.1, 63 = 0.08, 
bi — 0.05. Cutoff radius Rcut changes from 0.5 to 7 
with step 0.01. With increasing Rcut, the trajectories 
move to small Re(2), |Im(z)l (with the exception of 
the spiralling trajectory). 

Or 



mation functions only, and one fourth order trans- 
formation. The first order transformations keep un- 
changed the zero initial potential but affect the hyper- 
bolic transformation functions, producing only shifts 
in their arguments. Thus the potential is calculated by 
Eq. ^ with Vo = and the fourth order Wronskian 
W{vi, . ..,V4). Here Vi = sinh(6jr - 77^), r]i = Q + Q 
where for i = 3,4 are calculated by Eq. (^0]) and 
Q for i = 1 , . . . , 4 are calculated by the same formula 
with fli and 02 replaced by 03 and a^. After some 
calculations we obtain an explicit expression for the 
Wronskian, 



i—l J— 2+1 

x{bl - 6f ) cosh^i cosh^j sinh^fc sinh^;, (46) 

where — biV — rji and, in every term, k > I take the 
values from the set (1, 2, 3, 4) different from the values 
of i and j. An explicit expression for the obtained 
potential is rather involved and we omit it here. Its 
typical behavior is shown on Fig. [U solid line. 

The Jost function for the two resonance potential 
follows from Eq. (fTT 



F2res — Y\_ 



k + ibj 



(47) 



Thus the resonances occur aX k = ki = —02 + iai 
and k — k2 ^ —03 + ia^ with the mirror poles at 
k = a2 + iai and k — + ia^. 

The resonance behavior of the cross section is more 
visible when S'-matrix poles are close enough to the 



-10- 



-20- 



10 20 30 40 50 60 70 

Re(z) 



FIG. 9: Trajectories of ^-matrix poles (|35p for two- 
resonance potential at ai = —0.1, 02 = —2, 03 = 
-0.08, ai = -3, bi = 0.2, 62 = 0.14, 63 = 0.08, and 
&4 = 0.05. Cutoff radius Rcut changes from 0.5 to 7 
with step 0.01. With increasing Rcut, the trajectories 
move to small Re(2), |Im(2)| (with the exception of 
the spiralling trajectory). 



real axis. This is achieved by a proper choice of the 
Bargmann potential parameters. Results of our cal- 
culations are presented in Figs. [THll For the set of 
parameters chosen in these figures the complex ener- 
gies have the values Ei = (—02 -I- iai)^ = 3.99 — iO.4 
and E2 = (-04 + W3)2 ^ 8.9936 - iO.48. The figures 
show the same features as those obtained for the one- 
resonance case. Two narrow resonances stand sepa- 
rately from the chain of the cutoff poles. The poles 
of the S matrix are determined well enough when cal- 
culated according to Eq. p5|) . Solutions of the fixed- 
point approximation (j33p identify correctly the posi- 
tions of the physical poles if they are close enough to 
the real axis. 
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VI. CONCLUSION 

In the present paper, we considered the application 
of the FPO formahsm to potential scattering in the 
single-channel case. For this purpose we analytically 
constructed model potentials being a generalization of 
Bargmann potentials to resonance states with one and 
two resonances at given energies. 

In the FPO formalism, the corresponding spectro- 
scopic and scattering information is obtained from 
the non-Hermitian Hamiltonian Hqs and the S ma- 
trix derived by means of Hcs- The Hamiltonian Hcs 
describes the localized part of the system under the 
influence of the coupling to the continuum. In the 
present paper, it is obtained in the framework of the 
tight-binding model. 

First we compared the phase shifts obtained nu- 
merically in both methods. We received an astonish- 
ing good agreement of the results obtained in the two 
models. The phase shift notifies only the physical res- 
onances. 

To compare the spectroscopic values obtained in the 
two models, we are confronted with the problem that 
the eigenvalues of H^g involved in the S matrix are 
energy dependent while the exactly solvable poten- 
tials provide us only the poles of the S matrix. We 
therefore have to determine the poles of the S matrix 
also in the framework of the FPO formalism. The 
standard fixed-point approximation (j33p for the po- 
sitions of the resonances is inadequate for this pur- 
pose since the widths are determined by solving Eq. 
([M]) at the positions of the resonances. Hence the 



widths are, generally, erroneous. The results of our 
calculations show clearly that the fixed-point approx- 
imation gives, nevertheless, reasonable values if the 
widths are small enough. In order to determine ex- 
actly the poles of the S matrix in the framework of 
the FPO formalism, we solved the nonlinear equation 
dSlD, Det [E - Hcff {E)] = 0, in the complex E plane. 

The results for the spectroscopic values are the fol- 
lowing. The truncation of the potential in the tight- 
binding FPO model leads to the appearance of spuri- 
ous solutions of Eq. ([55]) just as in the well-known case 
of the S'- matrix cutoff poles j22|. The physical reso- 
nances of the truncated Bargmann potentials are well 
described by the complex energies satisfying Eq. (|^ , 
i.e., by the S matrix poles calculated in the FPO for- 
malism. Furthermore, the spurious solutions coincide 
with the cutoff poles of the scattering matrix. The last 
ones behave differently from the physical resonances 
in repeated calculations with different parameter sets. 
In our case the parameter is the cutoff radius Rcut 
(the potential is set to zero at coordinates r > Rcut)- 
The physical poles representing visual resonances are 
stable against variation of Rcut, in a certain range, in 
contrast to the cutoff poles that are not stable. 
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